Signatures of AGN model for UHECR 



V. Berezinsky 1 ), A. Gazizov 2 ' and S. Grigorieva 3 * 1 
x ' INFN, Laboratory Nazionali del Gran Sasso, 1-67010 Assergi (AQ), Italy 
B.I.Stepanov Institute of Physics of the National Academy of Sciences of Belarus, F.Skariny Ave. 68, 220062 

Minsk, Belarus 

3 ' Institute for Nuclear Research of Russian Academy of Sciences, 60th Anniversary of October Revolution 

Prospect 7A, 117312 Moscow, Russia 
(February 2, 2008) 

We demonstrate that the energy spectra of Ultra High Energy Cosmic rays (UHECR) 
as observed by AGASA, Fly's Eye, HiRes and Yakutsk detectors, have the imprints of 
UHE proton interaction with the CMB radiation in the form of the dip at E ~ 1 x 
10 19 eV, of the be ginning of the GZK cutoff, and of very good agreement with calculated 
spectrum shape. We argue that these data, combined with small-angle clustering and 
correlation with AGN (BL Lacs), point to the AGN model of UHECR origin at energies 
i?< lx 10 20 eV. Our consideration includes also the case when correlation with BL Lacs 
is excluded from the analysis. The excess of the events at E > 1 x 10 20 eV , which is 
observed by AGASA (but absent in HiRes data) can be explained by another component 
of UHECR, e.g. by UHECR from superheavy dark matter. 

PACS numbers: 98.70.Sa, 98.54.Cm, 95.85.S 



The nature of signal carriers and of the sources 
of UHECR are not yet established (for recent re- 
views sec fHH). The most natural primary parti- 
cles are extragalactic protons. Due to interaction 
with the CMB radiation the Ultra High Energy 
(UHE) protons from extragalactic sources are pre- 
dicted to have a sharp steepening of energy spec- 
trum, so called GZK cutoff jl. For uniformly dis- 
tributed sources, the GZK cutoff is characterized 
by energy Ei/ 2 where the integral spectrum cal- 
culated with energy losses taken into account be- 
comes twice lower than the power-law extrapola- 
tion from low energies, Eu2 = 5.3 x 10 19 eV ||. 

The particles with energies higher than 1 x 
10 20 eV are undoubtly observed. There are at least 
two "golden" events at energies 2 — 3 x 10 20 eV 
HQ] w ith very reliable energy determination (see 
also discussion in Ref. ||). However, the real con- 
tradiction with the existence of the GZK cutoff 
is observed only by AGASA (see Fig.l). Data of 
HiRes J§] are in good agreement with presence of 
the GZK cutoff (Fig. 1). The data of other detec- 
tors are not as conclusive, though a few events with 
energy higher than 1 x 10 20 eV are observed there 
(see [Q ) for a review and for recent discussion) . 

In this Letter we shall demonstrate that the ob- 
served spectra have the imprints of UHE proton 
interaction with the CMB radiation in the form of 
the dip at the energy E ~ 1 x 10 19 eV, produced by 
V + 7cmb — * P + e + + er interaction, in the form of 
the beginning of the GZK cutoff, and in the form 



of good agreement between predicted and observed 
spectra. We argue that at least at energies up to 
1 x 10 20 cV the data (spectrum, small-scale clus- 
tering [jn], and probably the correlation with BL 
Lac sources [jl2| can be explained in the model with 
AGN as the sources of UHE protons. 

Calculating the spectra, we shall use the cos- 
mological parameters as follows from recent ob- 
servations fl3|| : flat universe with f2 to t = 1 arid 
57a = 0.7. At small redshifts z, neutrinos, baryons 
and CDM behave as non-relativistic matter with 
Q m = 0.3. We shall use the Hubble constant 
Hq = 70 km/s Mpc. The relation between time 
and redshift is given by 



dt 



dz 



H (i + z)^n m (l + z)3 + n A , 
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The spectrum of UHE protons in the model with 
uniform distribution of the sources and with the 
power-law generation spectrum can be calculated 
using the formalism of Ref. IB, with the continuous 
energy losses from Refs. [ |14|Jl5| (note the difference 
in formulae due to cosmology with A term) : 
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where z m is a maximum redshift in the evolution 
of the sources, z g is a redshift at generation and 
E g (z g ) is energy of a proton at generation, for 
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present (z — 0) energy E; Cq = uqL p is cos- 
mic ray (CR) emissivity at z — (no and L p are 
space density of the sources and their CR lumi- 
nosity, respectively). As the general case we as- 
sume cosmological evolution of the sources given 
by C{z) = £o(l + z ) m , where the absence of evolu- 
tion corresponds to m = 0. All energies in Eq.(||) 
are given in GeV and luminosities in GeV/s. Dila- 
tion of energy interval is given by modified by 
Eq.(0) as 
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where bo(E) = dE/dt is the energy loss due to in- 
teraction with CMB photons at z — 0. Derivative 
db (E')/dE' at z = ( given in Ref. @) is taken 
at energy E' = (1 + z)E g {E, z). 

For particles with energies E > 1 x 10 17 eV, the 
maximum redshift for evolution of CR sources z rn 
is not important if it is larger than 4. In Fig.l the 
calculated spectra are compared with the data of 
AGASA, HiRes, Fly's Eye and Yakutsk detectors. 
Note that the theoretical spectra are the same in all 
four panels. We assume the generation spectrum 
cx E~ 2 at E < E c and cx E~^ at E > E c with a 
spectrum cutoff at E max . The calculations for this 
case are easily generalized from Eq. (Q) , as it is done 
by Eq.(10) from Ref. [Q. As parameters we have 
chosen £ max = 1 x 1CF eV, E c = 1 x 10 18 eV, 
7 9 = 2.7 for non-evolutionary (m = 0) model 
and 7 9 — 2.5 for evolutionary (m = 3) model. 
The required CR emissivity for both models is 
£0 ~ (2.5 - 3.5) x 10 46 ergs/Mpc 3 yr- The choice 
of E c is motivated by the observed chemical com- 
position. We assume that transition to extragalac- 
tic component occurs at the observed second knee 
g§- @, at E 2 « 4 x 10 17 eV. The observed rigid- 
ity cutoff for protons, E/Z = (4 — 5) X 10 15 eV 
according to KASCADE data corresponds 
to the cutoff energy for iron nuclei Em lxl0 17 eV. 
Remaining gap (1 — 4)xl0 17 eV can be filled by ul- 
traheavy nuclei with charge up to Z = 92 . The 
transition to the lighter chemical composition with 
a large fraction of protons has been observed at 
E > 3 x 10 17 eV by AGASA gl|, FE @, Yakutsk 
[|D and HiRes detectors. The data of HiRes 
[22 1 show that at E > 6 x 10 17 eV the protons can 
be the dominant component. 
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FIG. 1. The calculated spectra for non-evolutionary model (full lines), 
with parameters indicated in the text. Both curves were first normalized to 
adjusting the emissivity Co. To fit the data of HiRes, Fly's Eye and Yakuh 
factors 0.63, 0.80, and 1.7, respectively. 

Fig.l shows that that the signatures of interac- 
tion of UHE protons with the CMB radiation, the 
dip and the beginning of the GZK cutoff, are seen 
in the data. 

The most natural sources of the observed UHE 
protons are AGN. The required CR emissivity 
meet well the local emissivity of AGN, e.g. that 
of Seyfert galaxies is of order £sy ~ ns y Ls y ~ 
1 x 10 48 ergs Mpc 3 yr. The protons can be accel- 
erated in AGN up to energies of order ~ 10 21 eV 
p3[ |. An interesting possibility is acceleration in 
the jets by unipolar induction |24[ . The correla- 
tion of UHECR with BL Lacs i.e. with AGN 
whose jets are directed towards us, strongly sup- 
ports this mechanism. 

Two sets of observational data favor rectilinear 
propagation of UHE signal carriers in the universe. 

The first set is the above-mentioned correlation 
p2| with BL Lacs. The second one is the small- 
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angle clustering |TT|. Its most natural interpreta- 
tion is given |^5| in terms of statistically occasional 
arrival of two or three particles from a compact 
source. Such an interpretation needs rectilinear 
propagation of the primaries and large number of 
sources p5| . 

Since the proton origin of UHECR is almost 
proved, the correlation with BL Lacs directly 
implies the rectilinear propagation of UHE pro- 
tons. Thus, these correlations become supersensi- 
tive tools to measure extragalactic magnetic fields. 
Below we shall discuss what is the scale of this field 
and whether this field can be already excluded. 

Magnetic field must not produce the angular de- 
flection larger than angular resolution of sources in 
the detectors, which is typically i? res ~ 2.5°. The 
correlation is found in the energy range (4 — 8) x 
10 19 eV, for which the largest attenuation length is 
^att ~ 1000 Mpc. The required upper limit for the 
magnetic field, which is homogeneous on this scale, 
is B[ < 2 x 10~ 12 Z^) 00 G, where Ziooo is attenua- 
tion length for 4 x 10 19 eV protons in units of 1000 
Mpc. For a magnetic field with small homogeneity 
length Ihom the required upper limit is 

B < mcs ~ 6 x 10- 10 G, (4) 

ev^att'hoin 

where the numerical value is given for l at t ~ 
1000 Mpc and Z hom ~ 10 kpc. 

We argue that these fields are not excluded. 

The observed Faraday rotations give only the 
upper limits on large scale extragalactic magnetic 
field j26|. All known mechanisms of generation of 
the large scale cosmological magnetic field results 
in extremely weak magnetic field ~ 10~ 17 G or less 
(for a review see |27j]). The strong magnetic field 
can be generated in compact sources, presumably 
by dynamo mechanism, and spread away by the 
flow of the gas. These objects thus are surrounded 
by magnetic halos, where magnetic field can be es- 
timated or measured. The strong magnetic fields 
of order of lfiG are indeed observed in galaxies and 
their halos, in clusters of galaxies and in radiolobes 
of radiogalaxies. As an example one can consider 
our local surroundings. Milky Way belongs to the 
Local Group (LG) entering the Local Supercluster 
(LS). LG with a size ~ IMpc contains 40 dwarf 
galaxies, two giant spirals (M31 and Milky Way) 
and two intermediate size galaxies. The galac- 
tic winds cannot provide the appreciable magnetic 
field inside this structure. LS with a size of 10 - 30 
Mpc is a young system where dynamo mechanism 
cannot strengthen the primordial magnetic field. 



In fact LS is filled by galactic clouds submerged 
in the voids. The vast majority of the luminous 
galaxies reside in a small number of clouds: 98 % 
of all galaxies in 11 clouds pq ]. Thus, accepting 
the hypothesis of generation of magnetic fields in 
compact sources, one arrives at the perforated pic- 
ture of the universe, with strong magnetic fields 
in the compact objects and their halos (magnetic 
bubbles produced by galactic winds) and with ex- 
tremely weak magnetic fields outside. However, 
even in this picture there is a scattering of UHE 
protons off the magnetic bubbles and the scatter- 
ing length is ^ sc ~ l/irR 2 n : where R is the radius 
ofa magnetic bubble and n is their space density. 
Among different structures, the largest contribu- 
tion is given by galaxy clusters which can provide 
l sc ~ (1 - 2) x 10 3 Mpc. 

Leaving the correlation of UHECR with AGN 
(BL Lacs) as an open problem for future observa- 
tions we turn now to the alternative of excluding 
these correlations from analysis. The small-angle 
clustering then can be probably explained in the 
other extreme case of very strong magnetic field 
due to lensing effect ^9|. Rectilinear propagation 
of UHE protons is not needed anymore. 

Till now we discussed mostly UHECR at 
E<lx 10 20 eV. At higher energies, as AGASA 
data show, there might be an excess of events. 
They should be interpreted as the new compo- 
nent. One possibility is given by Superheavy Dark 
Matter (SDMP) @. The spectrum of UHE par- 
ticles, produced at the SHDM decay is now reli- 
ably calculated using the different methods pl| . 
All calculations give very similar spectrum with 
7 9 w 1.9 — 2.0. This spectrum is shown in Fig.l. It 
describes well the observed AGASA excess. The 
primary particles are predicted to be UHE pho- 
tons, which at energy E > 1 X 10 20 eV cannot be 
excluded by available experimental data, in par- 
ticular by the inclined Haverah Park showers |32j 

In conclusion, the observed energy spectra re- 
veal the signatures of interaction of UHE protons 
with the CMB in the form of the dip, beginning of 
the GZK cutoff and the good agreement with the 
predicted spectrum. Combined with small-angle 
clustering Q and correlation with BL Lacs | Q , 
these data require the rectilinear propagation of 
UHE protons and AGN as their sources. The cor- 
relation with AGN (BL Lacs) becomes thus the 
most sensitive method of measuring very weak 
extragalactic magnetic fields. In case this cor- 
relation is not confirmed by future observational 
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data, and the small-angle clustering is explained by 
some more sophisticated phenomena (e.g. by mag- 
netic lensing) , AGN remain the favorable UHECR 
sources, but as one of the possible candidates. 
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